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ABSTRACT: 



General formulas for the first-order effects of the 
oblateness of the Earth upon the orbits of satellites 
Were reported in 19S? by R. E. Roberson and D, G. 
King-Hele. Their derivations employed methods which 
are unfamiliar to those not having studied the mathe- 
matics of astronomy. In this Technical Report / 
Research Paper No. 76, an elementary derivation is 
given for these useful and important results. 



This paper will be presented 26 May 1967 at the 
Canadian Congress of Applied Mechanics , Quebec. 
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An Elementary Derivation of the Principal Effects of the 
Oblateness of the Earth upon the Orbits of Satellites 

by 

John E. Brock 

Professor of Mechanical Engineering 
Naval Postgraduate School 
Monterey , California 

Independently of each other, in 1957 R. E. Roberson in the 
United States and D, King-Hele in Great Britain provided general 
formulas giving the first order effects of the oblateness of the 
Earth upon satellite orbits, previous analyses by others having 
contained inconvenient restrictions. Roberson and King-Hele obtained 
their results by the use of methods which are familiar to those who 
are trained in the mathematics of Astronomy but not familiar to 
those having only the usual training in Engineering Mechanics, who, 
accordingly, do not have easy access to these useful and important 
results . 

In the present paper these results are derived by a quite ele- 
mentary procedure requiring nothing more than elementary geometry, 
vector algebra, and mechanics. An appendix considers a few simpli- 
fied models of the Earth which have the same first order effects as 
does an oblate Earth. This material should be of interest to en- 
gineers engaged in the space effort as a means of understanding 
the primary perturbations which result from oblateness. 

A complete table 6f notation is given in Appendix B. This 
table is not arranged alphabetically, however, since so many of the 
definitions of the quantities depend on others in the list. Instead 
the arrangement is in "operational sequence", that is, in the order 
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in which the symbols logically occur in the development. The 
notation agrees generally with that used in as tronautical litera- 
ture; for example, V means geopotential, V means true anomaly, and 
V denotes satellite (vector) velocity. 

First we establish some coordinate systems. 0 is the mass 
center of the Earth and n is a unit vector parallel to the axis of 
rotation of the Earth pointing North, a and b are fixed unit vec- 
tors in the Earth’s equatorial plane so as to form a right handed 
cartesian system a, b, n; actually, of course, because of the pre- 
cession of the equinoxes and other slow motions, these vectors are 
not truly fixed, but we may regard them as being fixed for the pur- 
poses of this analysis. We could take a to be in the direction of 
the projection upon the Earth’s equatorial plane of the line from 
0 to the first point of Aires, but there seems to be no uniformity 
in this matter. We observe that conventionally the orbits of Earth 
satellites are referred to the Earth’s equatorial plane in the same 
way that planetary orbits are referred to the plane of the ecliptic, 
but there does not seem to be perfect uniformity with regard to a 
reference position for longitude such as that provided by the first 
point of Aires. 

Let Q be the instantaneous position of the satellite, with r * 

OQ and V = r, the dot denoting differentiation with respect to time. 

We define the scalar coordinates of satellite position by the equations 
X = F*a, y * r*b, z * r*n [la,b,c] 

so that 

[ 2 ] 



r “ xa + i/b + an 
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The geopotential for an oblate Earth (see Appendix A) may be 
approximated by the expression 



where p » is the product of the universal constant of gravitation 
by the mass of the Earth, J is the coefficient of the second zonal 



satellite, the second term in the square brackets in Eq. [3] is much 
smaller than unity. For a further discussion of Eq. [3], see Appen- 
dix A. 

The force, per unit mass of satellite, exerted on the satellite 
by the gravitational attraction of the Earth is 



V * -pr ^(1 + Jr ^(1 - 3s^r ^)/3] 



[3] 



harmonic, and r = /r« r. An approximate numerical value for p is 

3 —2 

9.66x10^ (statute miles) '^(seconds) . At the surface of the Earth, 

-2 

at the equator where r = 3964 statute miles, the product Jr has 
-3 

the value 1.623x10 (dimensionless). Thus, evidently, for an Earth 




[4] 



2 2 2 

and to evaluate this expression, we note first that r = x + y + a 



so that 




9x dr^ 3x dr^ 3x 



/ /o\ / 2vn/2-l,^ . n-2 

= (n/2)(r ) (2x) = arr 



[5] 



with similar formulas in which y and a replace x. Thus 



f = ygradfr"^ + 




and by using Eq. [2] this may be more briefly written as 
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f - - yr - Jyr* ^r(l - - 2J\izr 

_2 

■ - yr • AJye^ - BJyn [6] 

where - r/r is a unit vector parallel to r and 

A » r»"^(l - , B = 2sr* ^ [7a, b] 

Note that if J were zero in Eq. [7] , we would have the simple inverse 
square law of attraction toward the mass center of the Earth. 

We now define h to be the angular momentum (with respect to 0) 
of a unit mass of satellite, i. e., 

h » rxv [8] 



Since 



V = f - d(re )/dt = re + re 

r r r 

and since e *e » 0, we have 
r r 

2 

rxh = rx(rxv) * r*vr - r*rv * r (e *ve -v) 

r r 

2 • • • 3 • 

= r [(r + 0)e - re - re ] « -r e [9] 

r r r r 

We now use the fundamental principle of Mechanics that the rate of 
change of angular momentum (about 0 which we consider to be a fixed 
point since we are concerned with the motion of Q relative to 0) is 
equal to the moment about 0 of the applied force; that is 

fi * rxf [10] 

In the present case with only one moving particle, this formula is 
very easy to obtain from Eq. [8] , from which 

h = fxv + rxv * rxv = rxf 



since if is the acceleration which is equal to f, the force per unit 
mass . 



Next, consider the expression Vxh. Taking the derivative with 
respect to time and using the notation ■ h/^, where h ■ h , we 
have 

d 

*j^(vxh) - vxh + vxh * fxh + vx(rxf) 
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_3 

= - pr rxh - Jp(Ae^ + Bn)xh + vx[rx(-DJpn) ] 

- pe^ + pj?z[Ae^xe^ + Be^xn + B?i~^vx(nxr)] [ll) 

Integrating this expression with respect to time gives 

vxh = p(e^ + s) [121 

where S is a vector whose time-derivative satisfies 

s = Jh[A6j^xe^ + BCj^xn + B;i“^vx(nxr)l [13] 

Also, from Eq. [10] and Eq. [6], 

h = BJynxr [14] 

Since we have seen, in Eq. [3], that the terras involving J are small 
compared to those not containing J, we regard J as a "small” quantity. 
Thus, Eq. [13] and Eq, [14] indicate that S and h are "nearly" con- 
stant vectors; that is, their magnitudes and orientations are but 
slowly changing. 

Next, consider the scalar quantity 

2 

r*vxh = rxvh = h ^ y(r +r*s) 

= yr(l + = yr»(l + s cosv) [15] 

where s = /| •$ and cos V ^ e^*S/s. V is the angle, shown in Fig. 1, 
which is known as the true anomaly. If we define the quantity 

p = [16] 

which has the dimensions of a length, we may write 

p =* r»(l + s cos v) [17] 

• — 1 * * 

In Eq. [17] only r and V vary "rapidly", p = 2h\i h and s both 
involving J as a factor. Eq. [15] is that of an ellipse, called the 
instantaneous elliptical orbit (lEO) , from which the actual orbit of 
the satellite begins to differ as time goes on. The point on this lEO 
(point P in Fig. 1) corresponding to y ■ 0 is called the instantaneous 
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perigee point, toward which the unit vector = S/s points. Also, 

the eccentricity, s, of the lEO is called the instantaneous eccentricity. 

The angle i (0 m i , satisfying cos i = is the (instantaneous) 

inclination of the orbit plane. A unit vector the direction ♦ 

i, e. , the vector j » being perpendicular to both the plane 

of the lEO and the equatorial plane, thus points toward the descending 

node D (Fig. 1) of the former; -j points toward the ascending node A. 

The angle o) * <AOP is called the argument of perigee. The quantity p 

is seen to be the semi-latus rectum of the lEO. 

We will use the unit vector i = such that i, j, e, form a 

n n 

right handed triad, slowly changing in orientation since j and are 
not constant. We let u and W be the angles BOQ and BOP, respectively 
(Fig. 1) so that cos u » i *6 , cos w = i *e . Note that u ^ V W and 

that W is ‘^nearly" constant since the lEO changes only slowly. This 

• • 

implies that u - V, approximately, a result which will be used repeat- 
edly. Note also that u = w + tt/ 2. We also introduce the unit vector 
that e^, also form a slowly moving unit triad. 

Finally, we introduce the unit vector I = jxn so that I , j , n is a unit 
triad which rotates slowly about n* Thus, in particular, we have 



i “ I cos £ + n sin i [18] 

i cos u + 3 sin u [19] 

A littiA manipulation easily establishes the following relations. 

nxr ■ (-1 sin u + j cos i cos u) [20] 

2 « re *n >• r sin i cos u [21] 

r 



Now it is obvious that, to the approximation with which we are working, 
we have 

e - [22] 

r hr 

but it is also easily possible to obtain this result formally as follows 
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• i 1 • • • 

* 1 cos u + j sin - (1 sin u - J cos u)u 



• • • • 

» (J cos u - 1 sin u)u 



since 1 and j nearly vanish. However, 

e, ■ e. sin M + e. x1 cos w » -i sin u + j cos u 
• h r h h 



and Eq. [22] is established. From this we have 




h 




[22a] 



[23] 



Substituting Eqs. [7b], [20], [21], and [23] into Eq. [lA] , expanding cos v 
■ cos(w - w) , and multiplying by At, we obtain 



It is elementary but tedious to obtain the indefinite integral of this 
expression and since the result of such an exercise does not seem to be 
of particular value, we will perform the integration for one complete 
passage only, as u increases by 2 tt,. to obtain what we designate as Ah, 
the net change in h per passage. To do this most easily, note that for 
non-negative integers M, N 



Since ‘j.ejj ■ 0, the variation, in the magnitude of h is of higher order 
in J. Thus we also have 



dh - 2Jpp ^ sin i (1 +s cos u cos w + s sin u sin w) 

2 

x(-I sin u cos u + j cos i cos M)dM 



[24] 




[25] 



unless both M and N are even. We easily obtain 

Ah * ^ sin 2i 



[26] 



-2 -1 -2 . 

A6j^ “ jjpir?: p sin 2i - jJiip sin 2i 

using Eq. [16], and from the definition of j we have 

Aj - A(6j^xn CSC i) ■ ACj^xn esc i + e^xn Aese i 



[ 27 ] 



I 
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-2 . . -2 . 

AJ • 2JiTp jxn cos i + jAcsc i - 2J7rp I cos i [28] 

since evidently we must omit the term in j because j*j • We may also 
see in another way that there is no j-component since the variation ini 
is of higher order in J, Ah being perpendicular to the plane in which 
angle ^ is measured. We observe that the line of nodes rotates (in a 
sense opposite to that of the satellite itself if i < tt/ 2) with mean 
angular velocity 

* -2pj7rp ^cos'^i [29] 

where P is the orbital period, which differs only by a term of order J 
from the value 



„ . 3/2 -1/2 

P » 1-na y 



[30] 



(a is the semimajor axis of ellipse) which would obtain if J were zero. 

We use the asterisk in this formula and later to indicate a mean value 
in the sense that we have indicated. Thus we can finally obtain the very 
important formula 

B* - -nc cos i [31] 



where 






[321 

-2 



in which rg represents the equatorial radius of the earth and ^ * vir^ 
represents the acceleration of gravity at the surface. We have used the 
following relations 

" ^(y-0) p/(l-e) - 2p/(l-8^) 

')T„«-2r,-l , -7/2,. 2,-2 1/2 T -7/2,, 2,-2 1/2^ -1 _ 

2J7Tp P « Ja (1-a ) y ■ Ja (1-a ) g v ■ C 

E 

Ratuming now to Eq. [13], we wish to express all vectors in terms of 
the (slowly moving) triad , e^, 6^. We have 



9 



e » e cos V + sin v: n * e, cos i + i sin i 

r s 0 " h 

1 * e cos ii?-e_sinL;; v=f=^ -f rue, xe 
s 0 ' r h r 

It is now a matter of straightforward, if tedious, substitution into 
Eq. [13]. Some intermediate evaluations are shown below. 



e, xe » 
h r 


-e sin U + 
s 


®e 


V 






e^xn = 

h 


e, xi sin i « 
n 


sin i 


(e sin w + 
s 




cos w) 


nxe = 
r 




e 

s 




"e 




% 






sin i cos 


w 


- sin i sin 


w 


cos i 






cos V 




sin V 




0 



= -e cos i sin U -f e. cos i cos i; -f 6, sin i sin u 
s 0 h 

• 2 • 

vx(nxr) = rre x(nxe ) + r w(e, xe )x(nxe ) 
r r h r r 

= rr(e sin i sin u sin v - e„ sin i sin u cos y + e. cos i) 
s 0 h 

+ r u sin i sin u (e^ cos y + sin v) 

2 • 

From Eq, [22a] we have r u ^ h and from Eq, [17] and Eq. [23] we get 

• 3 -1 • -1 

rr - r p su sin u = hs sin u (1 + a cos v) . Thus we can write 

Vx(nxr) = (1 + a cos v) ^^{sin i sin u[e^(a + cos v) + sin v] 

+ e, a sin V cos i] 
n 

-4 -2* 2 

We also note that we may write 3hr = Jp v{\ + a cos v) by using 

• • 

Eq. [23] and Eq. [17] and the fact that u * V approximately. Thus, upon 
substituting into Eq. [13], we find 

• -2 2.2 

S “ Jp (1+scos u){(l-5sin i cos y)(l + ecos u) (-6 sin V + 

s 

V 

2 , 

e„cos y) + 2(1 +s cos y) sin i cos u (e slnw+ e„cos w) + 

0 S 0 

2 ^ 

2 sin i sin u cos w [e (s+cos y) + e..sin y] + 

s 0 

86, sin 2t sin v cos w} v 

n 

Multiplying by it and expanding functions of u in terms of y and 
W, u) being regarded as constant, gives an expression for dS which 
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may be integrated for general V. However, it is the net change per period 

which is of greatest interest, and upon use of Eq. [25], we find that 

-2 2 . 

AS » Jsirp [(2-3 sin '^)Gq “ sin 2i sin u] [33] 

Since S*AS * 0, the variation of s is of higher order in J, and we have 

Ae^ * Jirp ^[(2- 3 sin^i)6Q - 6^ sin 2i sin u] [34] 

For the argument of perigee, we have 

- cos 0 ) * j *e * sin 
s 

- cos(arfAa)) * - cos o) + Ao) sin co 

» (j+AJ)-(e +A6 ) 

s s 

- j-e + j*Ae + Aj*e 

s s s 

=® sin W + Jirp ^(2-^3sin^i) j •Cq -f 2l*e^cos i] 

-2 2 . 

= sin w -h Jirp cos u (4 - 5 sin i) 

Since sin o) = cos u, we clearly have 

Ao) = jTrp"^(4 - 5 sln^t) [35] 

Dividing by P, the time of passage, we have 

u* = (C/2) (A - 5 sln^ i) [36] 

for the mean rate of change of argument of perigee. From Eq. [31] we see that 

•* 

0 * - C cos i 

is the mean rate of change of the longitude of the ascending node. Combining 
this result with that of Eq. [36] , we find that the mean rate of change of 
the longitude of perigee is 

w * (C/2) (4 - 5 sin^i - 2 cos i) [38] 

where to, the longitude of perigee, is the ’’angle'* to » + to measured in two 

different planes. 

The result given in Eq. [36] is sometimes incorrectly referred to as 
the average rate of rotation of the line of apsides; strictly speaking it is 
simply the average rate of change of argument of perigee. The average rate 
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of rotation of the line of apsides really is simply 



- ,(C/2)[(2 - 3 sin^t)6Q - 6^^ sin li sin w] 
the magnitude of which is 

B* - (C/2) [(2 - 3 sin^i)^ + (sin U sin 



[ 39 ] 



The difference between this result and the value given in Eq. [36] is known 
to astronomers as the difference between the draconitic and the inertial 
motion of the apsides (1). 

It is of some interest to have a formula for the mean angular velocity 
of the triad. It is not difficult to show that this is given by 

the expression 



B*riad = (C/2) [(2 - 3 sin^t)6j^ - 1 sin 2i] 



[40] 



The results of principal interest are given by Equations [26] through 
[40] . These equations or their equivalents seem to have been first obtained 
by R. E. Roberson (2) and by D. King-Hele (3). Both of them also discuss 
variations in the orbital period from the normal value 

P - 

but this analysis seems to be of a higher order of difficulty than that 
presented here. There are obvious practical difficulties in even defining 
an orbital period. It could, for example, refer to the successive passages 
through the perigee point (of the lEO determined at perigee) , or it could 
refer to successive passages of the descending node. The analysis leading 
to Eq. [12] has been adapted from a recent popular textbook (4). 

■' Grateful acknowledgment is made of helpful suggestions from 
Dr. W. E. Bleick of the Naval Postgraduate School and ffom Dr. W. W. 



King of Georgia Institute of Technology. 
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A large literature has been developed concerning the potential 
of various mass distributions. Much of it has been concerned with 
the problem of determining the potential at an exterior point of 
a homogeneous ellipsoid. This problem is surprisingly difficult, 
but if the figure is of revolution, the results can be expressed in 
elementary (even though quite complicated) terms . Among readily 
available treatments, that of Thomson and Tait (5) is quite readable 
even though it is quite old. In an editorial footnote added in 1912, 
H. Lamb refers to a history of the subject by Todhunter. The latter 
seems to have been reprinted recently, but the writer has not been 
able to find a copy (6). 

However, the Earth is certainly not a homogeneous ellipsoid, 
and, Indeed, the finer details of the mass distribution of the Earth 
are presently being revealed by sophisticated analysis of observations 
made on actual satellite orbits. Our purpose in this Appendix is 
merely to indicate the genesis of Eq. [3] and to present a very 
simple model which yields Eq. [3]. 

If 0 is the mass center of the mass distribution, Q the exterior 
point at which potential V is being evaluated, and T is the general 
location of mass element dm, we let r* denote the distnee OQ, p denote 
the distance OT, a denote the cosine of the angle QOT, and 6 denote 
the ratio p/r. It is well known then that the potential may be 
expressed in a series of integrals (over the mass distribution) of 
Legendre polynomials in a, viz.: 

V = -(G/r)/£P (a) 6" dm 
n 
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where the sum runs from zero to infinity and the symbol P( ) 

n 

denotes the Legendre polynomial of degree n. If the body is 
homogeneous and nearly spherical, with principal moments of 
inertia I . I . and I , this may be developed as 

V = (-M/r){l+(2r>^M)"hl +I +J -3(a^I +c^J )] + ...} 

X y z X y z 

where a, and q are the direction cosines of line 00 with 

respect to the principal axes (through 0). If we take the z 

axis as northward, and consider an oblate spherical Earth with 

semimajor axis A and semiminor axis B, this becomes 

2 2 

-Ir, . A -B ,, ^ 2 -2, . 

V = - yr» [1 + — (I - 3z r ) + ...] 

lOr 

and the coefficient J in Eq • [3] may be identified with the quantity 
2 2 

(A -B )/10. The deleted terms, indicated by ... in the two preceding 
equations indicate terms of higher order in r However, since the 
Earth is not a homogeneous spheroid, the value of J cannot be cal- 
culated from a knowledge of A = equatorial radius and B = polar 
radius, but, instead, must be inferred from observed perturbations 
of satellite orbits. The forgoing has given an explanation of tlie 
source of Eq . [3] without, however, offering a model of the Earth 
from which Eq . [3] may be derived. 

Another point of view is given in what follows. The oblateness 
of the Earth may be regarded as resulting in a deficiency of mass 
near the poles as compared to a homogeneous spherical Earth. Thus, 
we are led to consider a model composed of one positive point mass 
M located at 0 and two equal negative point masses (-/t?) located on 

the polar axis at equal distances d from 0; see Fig. 2. The potential 

-1 -1 ' -1 2 2 

V at Q is given by -V/G = Mr - mr^ - mr^ . Now r* = x + ^ so 

that x>7e have 

(r^ 2^^ = = r^±2zd+d^ = r^[l+r ^(d^^±2zd)] 
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where the upper (negative) sign refers to anc^ the lov/er (positive) 

sign refers to Combining these relations, we find 

Vr»/G^ = -(/.Vn) + 2 - (d/r)^(l - !>??' jv') + ... 

2 

plus terms involving higher powers of {d/r) . Since we propose to 
keep {d/r) quite small, we will he justified in neglecting them. Tlius 
writing n = a positive number, we have 

V = - - 2n) - n'/Gc?^r"^(l - + ... 

Now M{\-2n) - = M ^ the mass of the Earth, so that = V /(l-?n) 

h h 

Thus, recalling that p = GAf , we liave 

hi 

V = -pr^{l + [nd" / {\-2n) ]r ^'(1 - ^) + . . . } 

2 

so that we identify J/3 with the quantity nd. /(l-2n). Thus 
5.41x10'^ = = [n(l-2n)] 

and 

hi 

There is, of course, no unique decomposition. If we take n = 1/3, 

then m = M- and A/ = 3W with d = 0.023r_ = 92.3 statute miles, and 
EE’ E 

so on. The following tabulation gives a few such choices. 





M/M^ 


m/Vp. 


(d/r^)^ 


d. (miles) 




101 


50 


1.08*10“^ 


13.0 




21 


10 


5.41-10“^ 


29.2 




5 


2 


2.71*10~^ 


63.6 




3 


1 


5.41*10"^ 


92.3 




2 


.5 


1.08-10"^ 


130.4 


Since these values of 




are so small. 


then , a forteviovi. 


values of (d/r) 


will be even 


smaller, where 


r is the distance 


from the center 


of the 


earth 


to a satellite 


; thus, truncation 



- 1/2 

of the binomial series expansions of (r^ well justi- 

fied and any of these models does a good job of representing 



the Earth’s actual potential field. 
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Appendix B. 

List of Notations. 

Except for tlie division into general catagories (general, points, 
vectors, scalars), the listing is in the sequence that the symbols are 
introduced in the text rather than in alphabetical order. 

i4. General 

A denotes net change (per orbit passage) of quantity following. 

— denotes differentiation with respect to time, t. 
d t 

•k 

emphasizes mean value over one orbit passage. 

B. 'Points 

0 = mass center of the Earth. 

0 = position of satellite. 

P = perigee point of lEO. 

D = descending node of lEO. 

A = ascending node of lEO. 

B = point on orbit on vector 1 extended. 

T = location of mass element dm. 

C. Vectors 

n = unit North vector. 

a,b =: fixed unit vectors in Earth’s equatorial plane, 
r = vector 00 = position vector of satellite. 

V *s r = velocity of satellite 

^ = attractive force per unit mass of satellite. 

= Tfr = unit vector in direction of satellite position. 
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h = 


r><V =* angular momentum about 0 per unit mass of satellite. 




= \\lh = unit vector in direction of h. 


s = 


vector satisfying Eq. [13]. 


e 

s 


= S/s = unit vector from 0 toward P. 


j = 


e, xn CSC i = unit vector from 0 toward D. 
h 


i = 


jxe^ = unit vector of triad i, j, 




= e« xe = unit vector of triad e , G., 6, . 
h s s ’ 0 ’ h 


I = 


jxn = unit vector of triad I, j, n. 




= mean angular velocity of line of nodes. 



1 j, - angular velocity of e , 6., 6, triad. 

Triad s’ 0’ h 

Z). Scalars; dimensionality is indicated following semicolon. 



X = 


r*a; L. 


y = 


r*b; L. 


z = 


r*n; L, 


V = 


2 -2 

geopotential at point Q; L T . 


y = 


GM^ = 9 . 66x10 ^ (statute miles) (seconds) \ L T '. 


G = 


coefficient of universal gravitation; L^F 


J = 


coefficient of second zonal harmonic of potential 

-3 2 2 

1.623x10 L . 


r = 


/r-'r = distance OQ : L. 


A = 


r ^(1 - ^) ; AlsOy in Appendix Aj letters 


B = 


r i A and. B denote semi-axes of 

Izr ; L . ^ the Earth. 


t = 


time; T. 


k = 


/h *h = magnitude of angular momentum per unit mass; ^ . 


V = 


true anomaly, angle POQ = cos ^(e^-G ); dimensionless. 


P = 


2 -1 

h \i = semi-latus rectum of lEO; L* 


s = 


/s •$ = eccentricity of TEO; dimensionless. 
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i = inclination of plane of TEO with resnert to equatorial plane 

= cos dimensionless. 

U) = ar.(;ument of perigee = angle AOP ; dimensionless. 

u = angle BOO = cos ^(1-e^); dimensionless. 

w = angle BOP = cos ^(i*e^); dimensionless. 

a = semimajor axis of lEO; L. 

3/2 -1/2 

P = 2*na u = period of orbital passage; T. 

= coefficient in several important results: \ 

= radius of the Earth; L. 

-9 -2 

^ ^ ^ V ~ acceleration of gravity at surface of Earth; LT 

. * -1 

= mecm rate of change of argument of perigee; T 

•* -X 

Q = mean rate of change of longitude of ascending node; T 

^ -1 
= mean rate of cliange of longitude of perigee; T 

* “1 
- mean rate of rotation of line of apsides; T 

p = distance OT to mass element dm: L, 

a *= cosine of angle QOT ; dimensionless. 

B = o/y^y dimensionless. 

P^( ) = Legendre polynomial of degree n; dimensionless. 

2 

I .1 .1 = principal moments of inertia at 0; FLT . 

z ^ ^ 

a^hyC = direction cosines of OQ v/ith respect to principal axes; 
dimensionless . 



A,B = semiaxes of oblate Earth; L. 

2-1 

M = postttve mass at 0; FT L 

2 -1 

-m = negatzve mass on polar axis ; FT L 
d = distance from to -m: L. 
n = m/M\ dimensionless. 

2 -1 

M- = = mass of Earth; FT L 

E 
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Figure 1 Details of orbit of Earth satellite. Fieure la 
is looking south along axis of Earth : "Figure lb 
is looking along line of nodes from ascending node 
A toward descending node D; Figure Ic is true pro- 
jection of orbit. P is the perigee point, 0 is the 
center of the Earth, and Q is the location of the 
satellite. Angle AOR is a right angle. Q is the 
longitude of the ascending node, o) is the argument 
of perigee, V is the true anomaly, and angles w and to 
are as indicated. 
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Figure 2 Simplified model of the Farth 
consisting of positive mass 
at center and two negative masses 
on polar axis. 
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